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This paper describes a new method of defining the total aberrations of an optical system 
and its application to lens design. A system of partial differential equations defining optical 
image formation is written in the form of 4 by 4 matrices and is applied to the derivation 
of lens aberrations. It is shown that the form of the matrices facilitates the computation 
of these aberrations and that equations of condition on the lens parameters can be derived. 
The use of this method is illustrated through its application to a simple, well-known problem. 

1. Introduction 

In the course of designing a, lens, one invariably computes aberration functions whose 
magnitudes reflect qualitatively and quantitatively the performance of an optical system. 
Unfortunately these quantities are rather difficult to calculate, usually requiring a long program 
of ra}^ tracing, and, more importantly, are not readily correlated with the lens parameters. 
Usually a lens designer must rely on his experience when determining whether a change in 
a certain curvature, thickness, or index of refraction will affect a given aberration. 

In this paper a method of attack will be described which will at least mitigate these dif- 
ficulties. The application of this method results in the definition of a set of aberrations to- 
gether with a systematic scheme for their computation expressed as equations in which the 
system parameters appear. Yet the method is sufficiently flexible to meet the special require- 
ments of aiiv lens. For example, one system of aberrations can be defined for use with a copy 
lens and a second system can he constructed for application to an afocal system. 

Consider an optical system consisting of a system of lenses and fixed reference planes in 
object and image space which are not necessarily object and. image planes. Right-handed 
coordinate systems are chosen in object and image space, the coordinates of which are desig- 
nated x, ?/, :. and .r\ ?/', z' , respectively. We take the positive z and z' axes in the general 
direction of ray propagation and x-y plane and x'-y' plane coincident with the object and 
image reference planes. Let the reduced direction cosines of a ray in object and image space 
be £, 77, f and £', r/, f, where £ 2 + 2 r7 + f 2 =N 2 ., f ,2 +^ 2 +f /2 =A r/2 , AT and N' being the indices of 
refraction of object and image space. 

Associated with each ray passing through the optical system are two sets of four independ- 
ent coordinates, x, y, £, rj, and. x' , y- ', £', rj' '; for a given optical system either set completely 
determines a ray. In an obvious way a given optical system can be said to define x', y', £', and 
V as functions of x, y, £, and 77. In what follows we tacitly assume that in any optical system 
these four functions are continuous and possess continuous first-partial derivatives. 

Suppose every ray emanating from an object point converges to a single point on the image 
plane. Then x/ and ?/ ', considered as functions of the object space variables, arc unaffected by 
a change of £ and 77. Thus x f and 1/ are independent of £ and 77 and 

^= a; i=»« = »J=° (!) 

where the subscripts denote partial derivatives. The converse is also true; if there exists a pair 
of values for a; and y for which (1) is satisfied, then all rays emanating from the point (x, ?/) in 
object space converge to a single point (./;', y') in image space. Thus the four equations of (1) 
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provide a necessary and sufficient condition for a point to be imaged perfectly. By a similar 
argument it can be seen that in a region where (1) is satisfied a necessary and sufficient condition 
for a lens to be free from distortion is that 



tz=yl = ™>, Xy=yx=0, 



(2) 



where mis a constant and represents the magnification of the system. 

These two sets of equations provide conditions assuring perfect, distortion-free image 
formation. By the same token they may be taken as measures of the failure of an optical system 
to achieve perfection; that is to say, they are measures of the aberrations of the system. What 
follows is an elaboration of this basic idea. 

2. The Fundamental Optical Invariant 

The functional relationship described in the preceding section is a general transformation 
of which optical image formation is a special case. It is desirable therefore to impose restric- 
tions on the functions; to exclude from our attention those relationships that do not apply to 
optical systems. There are several ways in which this can be done. For the purposes of this 
study the most convenient method is through the use of Herzberger's fundamental optical 
invariant. This leads to a system of partial differential equations which must be satisfied by 
the functions x' , y f , £' , and -r\ if they are to describe the image formation of an optical system. 
Following Herzberger [1] these are written, 

Xx&+yzyy—x' y &—yyri' x ==0 
^+y f xV^—^x—y f ^x=^ 
x'xi'n +yWr,—K£'x —yW* =o 
*j£« +vW* — x> £y —yWv =° 
sjfi +yWr, — x '£v —yWv = 1 



(3) 



It is particularly convenient to express these in matrix form as 



where 



X= 



XJX=J, 

I x JCy Xj; Jiyj 

y'x y'v y\ y\ 

kx £y S| St? 

^ Vx Vy f]'t V 



(4) 



* 


— 1 





0^ 

1 





and A represents the transpose of X. 

It is important to note that X is a Jacobian matrix. An interesting corollary to (4) is that 



l*i=i, 



(5) 



or, in words, the Jacobian of the transformation representing optical image formation is equal 
to unity. 
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It is interesting to note that this corollary is equivalent to a theorem proven byStraubel [2] 
as long ago as 1902. He applied this to several problems including photometry. More recently 
Toraldo di Francia [3] used this theorem in describing a generalization of geometric optics. 

The X-matrix is particularly useful. First of all its elements consist of the partial deriva- 
tives of image space variables with respect to object space variables and therefore can be said to 
represent the geometrical aberrations of an optical system. Moreover, it appears as an indeter- 
minate in the equations derived from the fundamental optical invariant. Finally it can be de- 
composed into a product of X-matrices each representing a component of the lens. Thus, if we 
consider an optical system consisting of two lenses arranged so that the image plane of the 
first coincides with the object plane of the second then 



X=X\X\ 



(6) 



where X x is the A^-matrix of the first lens, X 2 of the second and X that of the composite op- 
tical system. 

The proof of (6) follows at once from the properties of Jacobian matrices and the rules for 
differentiating a function of a function. Let the first lens be defined by the four functions x, y, 
£, and rj of x ) y, £, and rj and let the second be defined by x', ?/, £', and 77' of x, y, f, and rj. The 
composite lens is then defined by x/ , y', £' , and rj' as functions of x,y, £, 77. Then differentiating 

x' x = ^Xx+y'-yx+x'-^ + xfi x 



Then 



x= 



= X 2 Jii. 



y'x y' v v\ y' v 



S/ %» £t £77 



X V £ 7? 

vi y'v v'i vk 



tj % «? K- 



r- - - - ■> 

x x Xy yi x 



ih Vu lit y, 

sx s.v s£ si? 

Kjfjx Vv Vz VvJ 



Clearly this can be extended to lenses consisting of any number of components. 

3. Transformation of the X-Matrix 



So far we have been dealing with the variables x, ?/, £, and f. It is convenient to apply 
a transformation, T, mapping these variables into a new system of coordinates, say a, ft 7, and 
5, given formally by writing the latter as continuously differentiate functions of the former. 
The same transformation is applied to both the object and image space variables. 

If we designate the M-matrix by 

ft Pfi Py Pi 

l'a lp 7y 7* 

vX 5'b 8' b' h J 



M= 



then 
where 



M=A'XA~ 



A= 



a x a u as, a v 
Pz ft ft ft 
lx 7y T€ Tt, 
A by ^ Si j 
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A r represents the same matrix but in image-space variables. Applying this transformation 
to equation (4) we have 

(A- 1 XA')A'- 1 JA'- 1 (A'XA- 1 )=A- 1 JA~ 1 , 
which reduces to 

M(A'- l JA'- 1 )M=(A'- 1 JA- 1 ). (8) 



If the transformation, T, is restricted by the equations 

AJA=J, A'JA' = J, 
MJM=J, 



then (8) becomes 



(9) 
(10) 



which is identical in form to (4) except that the variables are changed from x, y, £, and r\ to 
a, P, y, and 8. Thus we can change our variables in any way that is convenient, provided 
that (9) is satisfied, without adding to the complexity of the fundamental optical invariant. 

4. Rotational Symmetric Systems 

In the case of rotationally symmetric systems we are in a position to make a particularly 
advantageous choice of one of the new variables. In such a system the skewness invariant 
holds and 

p=x'7i'—y'£'=xri—y£. 

By letting y=xrj—y^=2), and y' = x'ri'—y'£', the M-matrix degenerates to 



M= 



r a a a'p a' p a s ^ 

Pa P'fi Pp Pi 
10 

Ja 8'fi K §'hJ 

When this is substituted into (10) we obtain the transformed fundamental optical invariant 

a a = 1 J rP p 8 a — fi a 8 v 

a^PX~P?p [ (11) 

a' 8 =P' p 8' s — P' s 8' p j 

P'J'p-P&=0 " 

P'J'i-fi&'a=0 V (12) 

P'/ 5 -P' s 8;=l. j 

We may consider the first two equations of (12) as a linear, homogeneous, simultaneous 
algebraic system in fi' a and 8' a . A necessary and sufficient condition for a non-trivial solution 
to exist is that the determinant of coefficients must vanish identically. However, the third 
equation of (12) show^s that this determinant is identical to a non-zero constant. We conclude 
that p' a =8 a =0. Therefore the first equation of (11) reduces to a a =l. The remaining equa- 
tions of the transformed fundamental optical invariant become 



a0=Ppdp — Ppf>p 
<X8=Ppdd—Psdp? 



(13) 
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and 

Solving (13) for fi' p and 5^, using (14), we obtain equivalent expressions 

fi p =P' B oL h —fta a 
8 p =8' B a 6 — 8' 6 a' B . 

Thus, under the assumption of rotational symmetry the Af-matrix reduces to 

"1 a'p a' p a b " 



M= 



^ P' p ft 

1 

^0 8^ 8 P 8' 8J 



By applying (15) this can be factored into 

C\ (H f\ a' a' «P 



M= 



o p; o ft 

10 



1 a h 

1 
— a* 1 



(14) 



(15) 



(16) 



5. Further Simplification 

Suppose we have a rotationalk -symmetric optical system defined by the functions a, 
$' and 5' of a, 0, p and 6, consisting of two components, the first defined by a, j3', and 6' of 
a, |3", p, and 5 and the second by a, (3, and 5 of a, 0, p and 6. Of course p and p are identical 
numerically but play different logical roles which must be kept distinct. As indicated in (6), 



"1 a' B a j, a b " 

fa is; # 

1 

^0 8'p 8 P 8 U 



1 a- a- a- * 
B v & 

% % % 

10 

,0 8'- 8- bj. 



1 a# a /; a 6 

ft ft, ft 

1 

.0 8b ~8 p ~hu 



(17) 



By interchanging the two right columns and the two bottom rows of each matrix in (17) we 
obtain 



"1 a' B a 8 a'^ C\ a- a- a p 

^ ft P' p 



8p 8' 5 8 P 
LO 1 



r B ^5 ^P 

8- a 8'- 8- 
1 



O a$ a s a p ~^ 

ft A /§, 

<5(-J ^ <5 7 , 

Ij) 1J 



(18) 



Now let 



A' = (a^), £'=(<4), 



with .4, #, r, /> and .1, 7^, C, D defined similarly. Then, by partitioning the matrices of 

(18) we obtain 
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1 


A' 


e 


C 





f 



B'^ 

D' 

1 



A 
C 

i 



D 

1 



A 
C 

f 



B 
D 

1 



(19) 



where e=( ft ) and /= (0,0). Then from (19) we can obtain 



C' = CC 



or 



ft 
h 



It has been shown that a consequence of (11) is that a'„=l, whence 

a'=a+\'(Ji,p,8). 
Applying this relation to the two components of the system under consideration, 

a ' = a+l(p,p,8) 
a =a+\(p,p,5). 



(20) 



(21) 



Thus 



-a=X' = X+X. 



(22) 



Thus the geometrical properties of a rotationally-symmetric optical system can be described 
completely by the system of equations 

(23) 



r=l 



and 



\8' 8 f J i-AtoJWi-i toi/tot-i) 



(24) 



where X * +1 is a function of /?*, £> and 5<, where p' = (} n , 8' = 8 n , and where P= Po, 5=5 . This 
is a system of five scalar equations of which four are independent by (14). 

If we replace the left members of these equations by functions representing a perfect 
optical system we are left with equations of condition on the parameter occurring in each of the 
components of the system. Alternatively a measure of the failure of the equations to be 
satisfied can be taken as a measure of the aberrations of the optical system. 

6. Functions for a Perfect Optical System. 

A possible choice for the variables, suitable for work with an optical system with object 
and image planes located a finite distance from the lens is 



a=arctan y/x^ 






(25) 



^fx Y +y 2 
It can be verified that these satisfy (9) when y=p. 

By integrating the equations for a perfect optical system (1) and (2), we obtain x' = rnx, 
y'=my, the constant of integration being eliminated by a suitable choice of coordinate axes. 
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From this 

a=a, and (3'=m/3< 

whence fi'p=m and /3 r 5=0. Substituting this in (14) we obtain 



Integrating: 

where ^ is an arbitrary function of 
variables in (25) we have 



Thus for a perfect optical system described by the 



-a=^\i = 0, 



and 






(26) 



(27) 



as our equations of condition. It can be shown that \p& when 0=0 equals the reciprocal of the 
focal length of the perfect lens. 

The complexity of the problem may be reduced somewhat by the following considerations. 
Suppose (27) were satisfied for meridian rays only. Then 

8'=8/m-\-6 2 (p,d,p) 

where €i=€ 2 =0 when / > = 0. If (26) is subsequently corrected for skew rays then a /3 = cx's = 0, 
which by (15) implies that p' p =8' p =0. Thus the derivatives of t\ and e 2 with respect to p 
are zero. Therefore these two functions are independent of p. Since t\ and e 2 vanish for a 
particular value of p } they must equal zero for all values and (27) is satisfied for skew rays as 
well as meridian rays. Thus, the solution of (27) for meridian rays and the solution of (26) 
for skew rays assures the correction of the optical system for skew rays. 

7, Transfer and Refraction Calculations 

The ultimate components of an optical system will be considered here to consist of (a) a 
transfer across a space between two planes and (b) a transfer from a plane (which will be called 
the reference plane) to a refracting surface, refraction across the surface and a transfer hack 
to the reference plane. The ray tracing formulas for the transfer are 

x'=x+t£/? 

y'=y+trj/t 



where t is the distance between the two planes, whence 



a'=a+arc tan 



tp 



=a-farc tanp//S5— arc tan />!ft'8', 

0' =>'(/?+ W/f) 2 +(^//Sf) 2 » 

s ,_ 8(0+ti/fi + (tplK)pl^ 
V(/W5/f) 2 +(^//3f) 2 ' 



(28) 
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For meridian rays the latter two formulas are 



J'=G8+«/f)l 



5' = 5, 



\ 



(29) 



J 



where f 2 =iV 2 — <5 2 . Accordingly the matrix representing transfer between two planes separated 
by a distance t } in a medium with index of refraction N, is 



m T - 



NHIf 







(30) 



The ray tracing formulas for a transfer from a plane to a refracting surface, refraction 
across the surface, and a transfer from the surface back to the plane are, respectively 

y=y+zvt, 
v' = ri+y4> 

y'=V-n'/!', . 

where ~z is a function of /3 denoting the refracting surface and where 



<t>-- 



;{ 



(r-5^)-V(r-5^) 2 -(i+^ 2 )( 



In terms of the a, j3, S coordinate system, these are 



N*-N' 2 )\- 



a=a-\- arctan 



zp 



]8=V(/S+2«/f >*+ (2>//3f )*' 

- 6(0+zS/t) + {p/fiGPlK) , 
V(/3+i5/f) 2 +(ijt>/^f) 2 '^ 



P'8' 



a'=a— arctan = — "^ — =a— arctan -^7+ arctan =-> 

/3 2 r'-s/35' 

/8' = V(is— 2«7f ') 2 + (ip/Pf ') 2 

s ,^j' (ii-i~8'/n-(p/0)(zp/pn . 

V(J-i5'/f') 2 +(i^r) 2 
38 



(31) 



(32) 



(33) 



(34) 



(35) 



Thus for skew rays 

n i) j) j) 

a'=a— arctan -^-7+ arctan 1^- — arctan ==+arctan — 

p pp(i- 8') , , V 

=a— arctan -^7+ arctan - — — — + a retail 7— 

=a— arctan -f-7— arctan = =— — = — [-arctan ™ ( :>('» 

^5 /3c5<M-(5 2 +^/0 2 ) ^ 

For meridian rays the formulas for £ and 5 reduce from those given in (33), (34), and (35) 



to 



+ m } (37) 



5 = 8 

i5'=j5 



"V (38) 



and 

where f 2 =iV 2 -5 2 and f' 2 =iV /2 -$' 2 . 



flv > w 



_ By differentiating the firsi equation of (37) with respect to /3 and 5, recalling thai I is a function 
of /3, we obtain 



so that 



where 



ft=(f/A)(^i/f 3 ), 



A=f— 27«. 



Th 



-=(f:f:)Kt wa) ( r /!4> ) 
=(? X *T) 

Dy differentiating* the equations of (38), 

where A' = f— 275'. From the definitions of A and A' and from (34) it can be shown that 
14-/ty(54-27f)/A'f=A£')/A'f, resulting in a simplification of m 6 . 



508&86- 59 4 
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Again, differentiating (39) we obtain 



m. 



\S'yS'v) \ 1 / 

J\ -zN'*lt*\ /A'/f 0\ 



When these three matrices are multiplied together, we have 
m R =m c m h m a 

(\ -zN' 2 l$ fZ \ /A'f/Af \ /l 2iV 2 /f 3 \ 

\° ! / W'z^/sjAA' AfVA'fj \0 1 /' 



(40) 



Note that the first and last factors are formally identical to m T and can therefore be combined 
conveniently with an adjoining transfer matrix. 

When the refracting surface is a sphere centered on the reference plane 

The minus sign assures that the sphere is concave to the right when c is positive, in accordance 
with the usual convention. Then 

Wb ~(o i Afr'^PAA' Af/A'fAo i / (41) 



where 



A=(0«+5f)/5, 

^ = 0§5 / +5fO/5 

0=— c 2 i(A— A'). 



A / = (05 / +ir)/^-'Y/A 2 -^-(^ 2 -iV /2 ); 



(42) 



8. A Simple Illustration 

As an example we consider the properties of a single refracting spherical surface. Con- 
sider an optical system consisting of an object plane, a spherical surface separating media of 
refractive indices N and N' and an image plane. Let the object and image planes be located 
distances t and t f from the center of the spherical surface. Figure 1 shows the case where t, 




/', and c are positive. From (27), (30), and (41) the equation of condition for meridian rays is 



/m \ 



1 N'H'/t"\/l -iA T,2 /f' 3 \/ A'f/Af 



1 /\0 



\/l zN>/i*\/l N*m*\ 
Af'/A'fA0 1 Ao 1 ) 

(43) 



For skew rays, by (26), (28), and (36), 

a' — a=— arctan 7^77+arctan —r— arc tan = L — — =-=0. 

0« ^ (384+(8 2 +p 2 /p 2 ) 

Expressing (43) in scalar form we obtain four equations. 

A'f/Af'+iV' 2 («'-2)^/c 2 i 2 AA / r' 2 =m ") 

tt' 4>/c 2 z 2 AA' = h 

N 2 (t+z)A' N' 2 {t'-z)A N 2 (t + z)N' 2 (t'-z)<l> 
Af 2 f + A'ff' 2 c^AA'f 2 ?' 2 

Af7A'f+A 2 (<+i)f>/c 2 PAA'f 2 =l/m.J 
We will pay particular attention to the third of these, which we write in the form 

N 2 (t + z)A'^' + N' 2 (t'-z)A^+(\/c 2 'z 2 )(t+z)N' 2 (t , -z)<j ) =0. 
Let 

(7 =Ai = /35+ir, 



(44) 



(45) 



(46) 



< T ' = A'-z = p8'+lf=j< r i -{l/c 2 )(N i -N'*). 



Also, 



where <£=— c 2 (<r— <r') by (42). Substituting these yields 

c 2 a^[N' 2 (t'-l)+N 2 (t+z)]-N^(t+z)(N 2 -N' 2 )+N 2 (t+J)lc i a 2 z-N i 'z+N'H']<l>=0. 
This may be written as 

c 2 <rH(N'H'+NH) + UN 2 -N' i )[c i <T 2 z-N 2 (t+I)]+N 2 (t+z)[c 2 a 2 z-N 2 z+N'H'}<t>=0. 



If 



Then 



N'H'+NH=0, 

[f(A/ 2 -A' 2 )+A 2 (<+I)<»][cVi-A^+i)]=0. 
We thus obtain two possible cases represented by the two equations 

UN 2 -N' i ) + N 2 (t+z)4>=() 
c i a i z-N i (t+z)=0. 



(47) 



(48) 
(49) 
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Whe»|8=5=0, 

1= -l/c,a= -N/c,<t'=-N'/c,{=N, and 4>=c(N-N'). 
Substituting in (48) we obtain 

N(N 2 -N' 2 ) + N 2 (t-l/c)c(N-N') = 0, 
N+N'+N(ct-1) = 0, 
t=—N'/Nc. 
Applying this to (47), 

t'=N/N'c. 



(50) 



(51) 



When /3 and 8 are not zero, it can be shown that (50) and (51) satisfy (48) for all values of 8 
provided that /3=0. Therefore only the axial point on the object plane is imaged perfectly. 
This result provides the aplanatic points of a sphere [9] (see fig. 2). 




Figure 2 



On making the same substitutions into (49) we find that /=0, which when applied to (47) 
yields 2' = 0. As in the previous case, these are true for all values of /3 and 8 provided that 
|8=0. This illustrates the well-known fact that the center of a sphere is imaged perfectly on 
itself [4]. 
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